HARMONIC ANALYSIS RELATED TO SCHRODINGER 

OPERATORS 



GESTUR OLAFSSON AND SHIJUN ZHENG 



Abstract. In this article we give an overview on some recent development of 
Littlewood-Paley theory for Schrodinger operators. We extend the Littlewood- 
Paley theory for special potentials considered in the authors' previous work. 
We elaborate our approach by considering potential in Cq° or Schwartz class in 
one dimension. In particular the low energy estimates are treated by establish- 
ing some new and refined asymptotics for the eigenfunctions and their Fourier 
transforms. We give maximal function characterization of the Besov spaces and 
Triebel-Lizorkin spaces associated with H. We then prove a spectral multiplier 
theorem on these spaces and derive Strichartz estimates for the wave equation 
with a potential. We also consider similar problem for the unbounded poten- 
tials in the Hermite and Laguerre cases, whose V = a\x\ 2 + b\x\~ 2 are known 
to be critical in the study of perturbation of nonlinear dispersive equations. 
This improves upon the previous results when we apply the upper Gaussian 
bound for the heat kernel and its gradient. 



1. Introduction 

The purpose of this article is to review recent development of harmonic analysis 
for differential operators, in purticular a Schrodinger operator H = — A + V, where 
V is a real- valued potential function on R". We are interested in developing the 
Littlewood-Paley theory for H in order to understand the associated function spaces 
and their roles in dispersive partial differential equations. 

This subject has been drawing increasing attention in the area of harmonic anal- 
ysis and PDE [2911311311311131131211122112313113311351123115!, to name only a 
few. The function space theory for H was introduced in [13 EU US] for the Hermite 
and Laguerre operators. In 01 [33J the authors considered Littlewood-Paley theory 
for H with special potentials in an effort to extend the function space theory to 
the bounded potential case. In this paper we will summerize and develop the fun- 
damental theory for general Schrodinger operators on R". Furthermore we obtain 
a Littlewood-Paley decomposition for L p spaces as well as Sobolev spaces using 
dyadic functions of H. We elaborate our approach by considering one dimensional 
H with V in <S(R), the Schwartz class. We will give outlines of the proofs for some 
of the main results and refer to the references, either old or new, for the detailed 
proofs. 
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1.1. Besov and Triebel-Lizorkin spaces. For a (Borel) measurable function cj>: 
K — » C we define by functional calculus 

/CO 
<K\)dE x 
-oo 

where H = J \dE\ is the spectral resolution of if. 
Let {^jj^oo C Cg°(R) be a dyadic system satisfying 

(i) supp tpj C {x : V~ 2 < \x\ < 2 j }, 

(ii) \ipf \x)\ < c k 2~ k i , Vj G Z, fc G No = {0} U N, 

OO 

(hi) J2 1^(^)1 ~ 1, V.t^0. 

j=-oo 

Let a € K, < p < oo,0 < q < oo and {ipj}j^z be as above. The homogeneous 
Besov space associated with H, denoted by B^ q {H). is defined to be the completion 
of <S(R n ) with respect to the quasi-norm 

oo 

(i-i) II/IIb- ( h) = ( E v aq \WAH)f\\D l/q ■ 

J=— oo 

Similarly, the homogeneous Triebel-Lizorkin space associated with H 7 denoted 
by Fp ,g (H), ft£M, 0<p<oo,0<<7<oois defined by the quasi-norm 

oo 

\\f\\^ HH) = \\( E v aq \^(H)f\ q ) 1/q \w- 

j=—oo 

We are mainly concerned with the following three interrelated problems. 

Problem 1.1. a. Littlewood-Paley theory for B(H) and F(H) using maxi- 
mal function characterization. 

b. Spectral multiplier theorem: Find sufficient condition for // € L°°(R) such 
that is bounded on L p , B(H), and 

c. Strichartz estimates for e ~ lt ^H and e~ ltH which measure the spacetime 
regularity of solutions to wave and Schrodingcr equations. 

The decay estimate in (|1.2[) has been known to be fundamental and useful in 
function space theory and spectral multiplier problem [T^l HU @1 [Ml GH] ■ We will 
see that it can be applied to characterize B{H) and F(H) spaces with full ragne of 
parameters < p, q < oo and show Mihlin-Hormander type multiplier result on L p , 
B(H) and F(H) spaces; for the multiplier problem we actually formulate a more 
general condition as in (|1.5|) . 

Let <j)(H)(x,y) denote the integral kernel of 4>{H). 

Assumption 1.2. Let <f>j E C^°(M.) satisfy the conditions in (i), (ii). Assume that 
for £ = 0, 1 and for every N G No there exists a constant cjy > such that for all 

2 (n+£)]/2 

(1.2) \V%{H){x,y)\ < ey ^— ^_ — . 

We will outline the proof of the fact in Section [4] that on the real line H satisfies 
Assumption 11.21 for V G C 3O (R). We discover that V being compactly supported 
and H having no resonance at zero are necessarily and sufficient for the gradient 
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estimate (£ = 1) in (|1.2p to hold in low energy — oo < j < 0. The kernel decay for 
£ = 0, j = was an open question in [38l B.7]. Q 

Define the Peetre type maximal function for H as: for j G Z, s > 

^3,sf{x) = SUp , — . 

J teR" (1 + W 2 \x - i|) s 

The following theorem gives a maximal function characterization of the homoge- 
neous spaces. 

Theorem 1.3. Suppose H satisfies Assumption \l.S\ and {fj} is a system satisfying 
(i)-(iii). The following statements hold. 

a) If < p < oo, < q < oo, a£l and s > n/p, then 

Wfh^iH) ~ \\{^U H )f}\U^) ■ 

b) If < p < oo, < q < oo, a£l and s > n/ min(p, q), then 

Wfh^iH) « l!{2 J >* s ( J ff)/}|| L!) ^ ) . 

It is well-known that such a characterization implies that any two dyadic systems 
satisfying (i)-(iii) give rise to equivalent norms on Bp q (H) and F^ q {H). Another 
consequence is that H has the lifting property. 



Corollary 1.4. Suppose H satisfies Assumvtion M.SX Let s,a G K and < p, q < 

oo. Then H s maps Bp' q (H) isomorphically and continuously onto B^~ s ' q {H). 
Moreover, \\H S f\\ga- s , q ^ H ^ ps || /II ^"''(/f) - The analogous statement holds for Fp' q (H). 

The proofs of Theorem 11.31 and Corollary 11.41 arc quite standard and can be 
found in [551E3]; see also [HI [15]. 

Using Calderon-Zygmund decomposition and Assumption 11.21 we can show that 
L p (R n ) = F^ 2 (H) if 1 < p < oo and obtain the Littlewood-Paley inequality for 
L p spaces. If in addition V G S, we can using lifting property of H to characterize 
the Sobolev spaces W 2a (R n ) = F^ 2 (H), the inhomogeneous versions of F^ 2 (H), 
with equivalent norms [55, 33;. 



Theorem 1.5. Let 1 < p < oo. The following statements hold. 

a) If H satisfies Assumvtion \1.2l then 

oo 

ll/IU-(E-)«||( E \Vj(H)f(.)\ 2 ) 1/2 \\ LHRn) . 

j=—oo 

b) If, in addition to the condition in a), V G 5(R"), then for all a G K 

oo 

\\f\\ w ?><R») « \mH).f\\L P{ ^) + ii(E 22ja i^( ff )/(-)i 2 ) 1/2 n^(«") > 

where {<5>, fj^^i is an inhomogeneous system satisfying (i), (ii) and (Hi 1 ) below. 



^For the sake of exposition we are not trying to pursue how singular V G L 1 (R) can be. Indeed, 
for each N the proof of Theorem l4Tl shows that (TO) holds for j > if \V^(t)\ < (i}"^- 2 " 6 , 
< s < N and that 1 11. 21 1 holds for j < if V € L 1 with compact support and H has no resonance 
at zero. 
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The analogous results above also hold for the inhomogenous spaces Bp' q (H), 
F£< q (H) if using the system {$, ^>j}f =1 C C£°(R) with supp$ C [-1, 1], Lp 3 satis- 
fying (i), (ii) and instead of (iii) 

oo 

(iii') |S(a:)| + 5^(1)1 « 1, Vs. 
i=i 

However, the homogeneous space, which contains both high and low energy anal- 
ysis of H, are essential and more useful in proving Strichartz inequality for wave 
equations; see e.g., (351 [30] and Section^ 

1.2. Spectral multipliers. For Problem 11.11 b. Mihlin-Hormander type spectral 
multipliers for H have been considered in [5^1 HPJ [TBI EES] arid more recently 3G . 

In the classical case Hq = —A, one can use Calderon-Zygmund lemma to prove 
the L p Fourier multiplier theorem by showing that the kernel of /i(-ffo) verifies the 
Hormander condition 

(1.3) / \K(x,y)-K(x,y)\dx<C 

J\x-y\>2\y-y\ 

if n satisfies certain smoothness condition; see e.g. [551131]. However, for a general 
elliptic operator (|1.3p is not available. For the Schrodinger operator H with V > 0, 
Hcbisch [5S] used heat kernel estimates (h.k.e) to prove a spectral multiplier theorem 
on L p . Later on the heat kernel approach has been further developed to deal with 
positive selfadjoint differential operators [13]. 

The question remains if the negative part of V is nonzero, in which case the 
upper Gaussian bound for e~ tH may not be valid. In [56] we are able to treat 
general V by replacing the h.k.e. with a (much) weaker condition for the pointwise 
decay of a spectral kernel, namely (|1.5[) . 

As in [56| the following hypothesis on H is the main ingredients in proving 
spectral multiplier theorem. Let 4>j(x) = (f>(2~~ 3 x). 

Assumption 1.6. a. (Weighted L 2 estimate) There exists s > n/2 such that 

(1.4) S uv\\\x-y\ s ct> j (H)(x,y)\\ Ll <c2W 2 - s W 2 Vj G Z, 

y 

where c = c(\\4>\\x^)- 

b. (Weighted L°° estimate) There exists a finite measure £ such that for all j E Z 

(1.5) \^{H){x, y )\<d [ y n i 2 {i + y/ 2 \x-y-u\y n -tdau), 

where d = c' (\\<f)\\ w n+e) . 

Here A* = {/ E A S (R») : ||/||x| = : su Pt>o ll/(^) r ?ll^' s < 00 }i where 77 is a fixed 
function in with support away from 0, R* = R\ {0} and X s is either W/(R) 
or C S (R), the Holder class [3PJHFJ. 

In one dimension Assumption 11.61 is satisfied for X s = C 1 if V E L},(M.) or 
V E L\(R) and H has no resonance at zero, where := {V : (1 + )' T V" E L 1 }. 
In three dimensions Assumption 1 1 . 61 is true for X s = ifV = V+—V- is in the 

2 Again we are not trying to seek optimal condition on V; one can show that if \d^V(x)\ < c^, 
|fc| < 2mo — 2 for some mo £ N, then b) is true for \a\ < mo. 
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Kato class with small Kato norm, namly ||V±||if < 2ir and H has no resonance at 
zero [57119]. 

Under Assumption 11.61 for H and the condition /i G X% for some s > n/2, the 
boundedness of n(H) on Bp q (H), 1 < p < oo is an immediate consequence of the 
LP result in [55]. To prove that fi(H) is bounded on F™' q (H) we use an LP(^) 
vector-valued version of the proof of the L p result by applying Calderon-Zygmund 
decomposition and the dyadic estimates (| 1 .4|) and (jl.5p . 



Ilfll^ 




ll/ll_p«.s (iJ ) : 






which 



Theorem 1.7. [33] Suppose H satisfies Assumvtion \1.6\ Let X s = C s orW^ and 
r] be a fixed function in (R* ) . // there exists some s > l| so that 

sup \\fi(t-)n\\ x ° < oo, 
t>o 

then n(H) is bounded on B£' q (H) and F^ q {H), 1 < p, q < oo, a € R. 

It is easy to observe that X = C s corresponds to the usual Mihlin condition and 
X = W% the Hormander condition. That the the exponent ^ is sharp has been 
noted in e.g., [41} [13]. Note that under the conditions in Assumption II. 61 which is 
an alternative condition than Assumption 11.21 Theorem 11.71 implies 

"(H) 

given any two system {^}j e 
as we have mentioned. Moreover, by interpolation and duality we obtain from the 
proof of Theorem 11.71 that 

F°- 2 {H)=L p 7 Kp<oo 

(see [33] for the inhomogeneous case) , which is part a) of Theorem 11.51 while under 
somehow more general conditions in Assumption 11.61 

Remark on Assumption 11.61 Assumption 11.61 is intrinsic in the sense that it 
only relies on the property of H and is independent of the multiplier /i. As can be 
seen from the proofs in [56] [34] , Inequalities (|1.4j) and (| 1 . 5[) are to control higher 
and lower energy estimates of fj,(H ) respectively. If in (|1.4[) letting £ = 5, the Dirac 
measure, then we obtain the following pointwise decay 

\^(H)(x, y)\ < c n , e 2^ 2 (l + V' 2 \x - 

which is only valid, in general, for nonnegative potentials. This is the reason why 
we call (b) a "weighted" pointwise estimate. 

The remaining of the paper is organized as follows. In Section ® we apply the 
interpolation properties of B(H), F(H) to obtain Strichartz estimates for e~ lt ^ . 
In Section [3J we provide the outlines of the proofs of Theorem 11.31 and Theorem 1 1.51 
under Assumption 11.21 In Section [4] we show that for any V in iS(R.), Assumption 
11.21 is verified for high energy and for any V in (R) , the low energy estimates 
holds in the absence of resonance. The proofs are based on certain new and refined 
estimates for the modified Jost functions and its Fourier transforms whose details 
are quite lengthy and will appear elsewhere. For unbounded potentials wc consider 
in Section[5]thc analogue of Thcorcm ll.3l for Hcrmite and Lagucrre operators, where 
V = a\x\ 2 + b\x\~ 2 , by using gradient estimates for e~ tH . Further, we would like 
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to mention that the literature in the area suggests that it is possible to consider 
analogous problems for H with (degenerate) magnetic potentials, cf. [481 EH El] 



1 " 



2 

3=1 



icijXj) , 



by following a similar approach developed here. 

2. Strichartz estimates for H 

It is well known that Strichartz estimates have useful applications in wellposed- 
ness problem for nonlinear dispersive equations [301 142[ |37j . Consider the following 
perturbed wave equation with a potential on 



(2.1) 



\ u tt + Hu = F(t, x) 
y u(0,x) = u (x), u t {0,x) = ui(x) 

whose solution is given by 

, , . r=. smttVH) /"* sin((i - s)VH) , , , 
u{t, x) = cos(tVH)u + — v ._ V + / — - ' — -F(s, -)ds. 

vH Jo vH 

When n > 2, the original Strichartz estimate for (|2.ip with V = reads [43] 
(2-2) llull 2n±l < c||/|| 

if i*o = 0, Ui = / and F = 0, where Wf : = W^iM™) denotes the homogeneous 
Sobolev space. We are interested in proving Strichartz estimates for (|2.ip in the 
Besov space scale, assuming that the dispersive estimate (|2.3j) is true. In this section 
for convenience we will use Bp' q {^/~H) instead of Bp' q (H), whose norm is given by 
(HID with H replaced by s/H. Note that B^'^VH) = B* q {H). 

Assumption 2.1. Assume that u(t,x), the solution to \2.1\) with uq = 0, u\ = f 
and F = 0, satisfies for all t =/= 

(2-3) ll«(*,-)||o <c|t|-("- 1 V3||/||. a ^ tI 

B 1 2 (VH) 

Dispersive estimates in ()2.3j) were obtained, for instance, in [3l for V being 
smooth and in [9j for V in the Kato class. 

The idea to treat (|2.ip is to combine the arguments in [3D] and [2S] for a free 
wave equation. The Littlewood-Palcy decomposition seems efficient in dealing with 
this type of estimates although we do not have available the scaling invariance for 
H or ipj(H), as is important and crucial in the classical case. 

For this purpose we will need some interpolation, duality and embedding proper- 
ties for Bp' q (^/ll), which are analogues of B^ ' q (M. n ) and can be derived as corollaries 
of Theorem 11.71 [531 HP] . Directly using the classical Besov spaces would encounter 
commuting problem. 

From the expression of u(t, x) we see that it is essential to estimate e ±lt ^~^ f . 
Note that (|2.3j) is equivalent to 



(2.4) oiloo^cltl-^- 1 )/^"- 1 )/ 2 !!^^)/!!!, 



where %(•) = ipj(vH)u(t, ■). Using (|2.4p and TT* argument in [25j[30] we obtained 
in [57] the following theorem. 
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Theorem 2.2. Let n > 2, q, r £ [2, oo] and s€l. Suppose 77 satisfies the estimate 

in 12.3]) . Then the following estimates hold. 

a) 

(2-5) W^fW^Bf'iVS) < 4f\\ B ^(VH) > 

where q — 2 and cr = cr(q, r) verifies the gap condition | + ^ = ^ — cr. 
&j Lei 7 C K &e an interval. 
ft 



. f sin((t - s)V77) , , , „ 
II J Q ^ J? ( s >-)rfs|| i? B ; ^ ( VH) < c \\ F hf B °+*°-^ {V Hy 

where q' , r' denote the usual Holder conjugate exponents of q, r, (q, r) ^ (2, 2 ™Z$ ) 
and (q,r) are wave-admissible, that is, - + < ■ 

The estimates in Theorem l2.2l are analogous to those in the case of zero potential 
[25] . Under additional condition, e.g., assuming 77 satisfies (|1.2p . we can release the 
restriction q = 2 in (|2.5p by applying Besov embedding inequality. Observe that 
then the Besov space method yields a sharper estimate than (|2.2p if in (|2.5p taking 
q = r = (2n + 2)/(n — 1), a = 1/2, s = and substituting 77^2/ for /, noticing 
that B° 2 {VH) ^ F?' 2 (y/H) = U provided r > 2. 

Problem: Does the endpoint estimate in Theorem 12.21 hold with (q, r) = (2, 2 ™Z% ) 
for n > 4? 

The endpoint estimates involving L q t 7^-norm were proved in |30| in dimensions 
> 4 (In dimensions 2 and 3 the endpoint estimates fail). We do not know the 
answer to the question in the problem for L q Bp 2 . Such a result would not only be 
sharper but technically might involve bilinear estimates that provide insight and 
deeper understanding of the non-scaling- invariant case when V ^ 0. We can also 
formulate similar result and problem for the Schrodinger equation with a potential; 
here we would rather refer to [231 [3TJ [331 [57] for further discussions. 

3. Outline of Proofs 

3.1. Proof of Theorem 11.31 As in [33] [19] or [49], Theorem 1 1.31 is a consequence 
of Peetre type maximal equality (Lemma 13. 2p and the well-known L p (^ 9 )-valued 
Feffcrman-Stcin maximal inequality The proof of Lemma 13.21 is standard and 
follows from Bernstein type inequality (Xemma l3.ip . Let 

™ JK ' tf £ (i + y/ 2 \ x -t\y 

Lemma 3.1. For s > 0, there exists a constant c n _ s > such that for all jeZ 

< c n ^J/ 2 ^J(x), V/ G 5(R"). 

Similar to [331 SS] Lemma 1370 can be easily proved using (|1.2p with N > n + s. 
Let M denote the Hardy-Littlewood maximal function 



(3.1) Mf(x) 



sup ii?T / \f(y)\ d y 

B3x JB 



where the supreme is taken over all balls B in M. 71 centered at x. 
Lemma 3.2. Let < r < oo and s = n/r. Then for all j E Z 
(3.2) ^ j J{x)<c n>r [M{W j {H)f\ r )] l,r {x), V/e5(I 
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Thus we have seen that the proof of Theorem 11.31 relies on the decay estimates 
in Assumption 11.21 In Section |4] and Section [5] we will prove such estimates for 
H = -A + V with V G S on M, and V = \x\ 2 and V = \x\ 2 + b\x\~ 2 , b > on 
R n , where for the latter, the Laguerre operator, initially defined on R", can be 
regarded as an operator acting on L 2 (M. n ) for "even" functions. 

3.2. Proof of Theorem Ol The same proof of the identification F°' 2 (H) = LP 
in the inhomogcneous case [331 Theorem 5.1] gives us 

(3.3) \\f\\F° :2 (H) ~ Wfhp , l<p<oo 

for any system {tpj}j£Z satisfying (i), (ii), (iii), by applying L p (^ 2 )-valued Calderon- 
Zygmund decomposition. □ 
For p = 1, Dziubahski and Zienkicwicz recently obtained a characterization of 
Hardy space associated with H using the heat operator. 

Theorem 3.3. ([T7]) Let V G L% +£ (W n ), n > 3, with V > being compactly 
supported. Then 

\\f\\n{, ~ H/H«i tom ~ lk/lli?i(R») - 

where = {/ € L 1 : sup t>0 |e _tff /(-)| G L 1 } and the weight iu is defined by 
w(x) = Hmt^oo J R „ e~ tH (x, y)dy. The norm in the atomic decomposition is defined 
as 

where the infimum runs over all representations / = J2j ^j a j> a j being Hy atoms 
satisfying (i) suppa C B(x ,r) := {x:\x- x \ < r}, (ii) Halloo < \B(x , r)| _1 , (iii) 
J a(x)w(x)dx = 0. 

Remark. It would be interesting to obtain a norm characterization for TL V of 
Littlewood-Paley type in the sense of Theorem 11.51 Note that if V ^ 0, then 
7i v ^ H 1 {W l ) 1 n > 3. Also, in the ID case and unbounded potential case (e.g. V 
being a (nonnegativc) polynomial), the p-atoms (p = 1) for Tt v (V ^ 0) are not 
variant of the local atoms, cf. pTTt [T5] . 

4. H=-£s + V, V in S(R) 

Associated to H there exists a decomposition L 2 (M) = H a c © 7~Lpp, where 7i ac 
is the absolute continuous subspace and TC PP the pure point subspace of L 2 (M). 
Let E ac , E pp be the corresponding orthogonal projections. If o- ac (H) denotes the 
absolute continuous spectrum and o~ pp (H) the pure point spectrum of H, then 
o~ac(H) — [0, oo ) and a pp (H) — {—A 2 } is a finite set of eigenvalues of H in (— oo, 0). 

4.1. Decay estimates of </>j(H). For A G R let e(x, A) = (1 + i?. (A 2 +iO)V)- 1 e lXx 
be the Lippman-Schwinger scattering eigenfunction and efe(x) the L 2 eigenfunction 
of H with eigenvalue —A 2 , where Rq(z) = (Hq — z) _1 is the resolvent of H = 
—d 2 jdx 2 with z G C \ [0, oo). Then if G Co, a continuous function with compact 
support, we have 

<KH)f(x)= [ K(x,y)f(y)dy 



HARMONIC ANALYSIS RELATED TO SCHRODINGER OPERATORS 



where K = K ac + K pp , 
(4.1) K ac (x, y) = (27T)" 1 / 4>{X 2 )e(x, X)e(y, X)dX 



is the kernel of 4>(H ac ), H ac = HE ac and K pp (x,y) = J2k ^(~^l) e k( x ) e k{y) is the 
kernel of <p(H)E pp . Let m, 112 be two linear independent solutions of (|4.3[) . H is 
said to have resonance at zero provided that the Wronskian of u\ , U2 vanishes at 
zero. 

Theorem 4.1. Let {fa} C C§° (R) satisfy 

a) If V G S , then for each I > and iV > there exists a constant cn,i so that for 
all j > 

(4.2) \difa(H)(x,y)\<c N , e 2^/ 2 (l + 2^ 2 \x-y\)- N . 

b) If V € C§° and H has no resonance at zero, then holds for each t = 0, 1, 
N > and all -00 < j < . 

Denote by Kj the kernel of 4>j{H) an d write i£j = Kj ac + Kj iPp . Since <7 PP is 
finite and according to e.g. [351 Theorem C.3.4], eigenfunctions of H belonging to 
n™ JV=0 (x)- N W 2 2m (R) = 5(R), ATpp satisfies g2J) trivially. Hence it is sufficient to 
deal with Kj >ac . 

Let R v (zj = (H- z)- 1 . Let W(X) be the Wronskian of /+, /_, then for A ^ 

( f+(x,±X)U(y,±X) 

Rv(\*±iO){x,y) = \ Mv Z%% t±x) J 

where f±(x, z) are the Jost functions that solve for > 

(4.3) - fi{x,z) + V{x)f±{x,z) = z 2 f ± (x,z) 



and satisfy 



j e ±lzx x — ► ±00 

J±{x,Z) — ► < x iia, TjKt 
777T C T C 



where t(,z),r±(z) are called the transmission and reflection coefficients [T2"] . 

Let m±(x,z) = e T%zx f±(x, z) be the modified Jost function. We obtain, from 
the resolvent formula of the spectral measure of H ac , that if x > y 

1 f°° 

(4.4) <t>(H ac )(x,y) = — J <j>{X 2 )m + {x,X)m-{y,X)t{X)e iX ^dX 

where t(X) = — 2iX/W(X); see e.g. [23]. One can prove that the above formula and 
(|4.1|) coincide. As a result the restriction x > y can be dropped. 

To analyze the kernel we need to study the analytic asymptotics of m±(x,k), 
namely the higher order mixed differentiability on m± concerning the smoothness 
and decay in x and k, the space and spectral parameters. The asymptotics of m± 
and their derivatives of first order were originally studied in [12] . However our 
estimates and certain formulas, including those of t(k) and the Fourier transforms 
of m±, arc more refined and delicate. 
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4.2. Derivatives of m(x,k). 

Lemma 4.2. Let V G S and £,n>0. Then m G C°°(M x K) and /or aZZ a, fc £ 

'(1 + max(0,Ta;)) n ^ > 1, 
(1 + max(0,Ta;)) n+1 ^ = 0. 



\d l x dkm±(x, k)\ < Cn t £ ■ 
lf£>0,n>l or I > 1, n > 0, then 



i fl i fl „ , M , . (1 + max(0, Tx)) r 



\k\ n 

The proof of the lemma are based on the integral equation [T2"] 

/>oo 

(4.5) m + (x,k) = l + h(t-x,k)V(t)m + (t,k)dt, 

J X 

where h(t — x, k) = J Q x e 2lku du, and the equation for its mixed partial derivatives: 
for 1= 1,2,3,..., n>0 

(4.6) did£m + (x, fc) = - E ■ W / e 2ik(f -*Kt ~ ^ ^(Vm^dt, 

j_Q \3 J Jx 

where rriAj = d l k m + (t, k). 

The first estimate in the lemma is also a consequence of Lemma 14.51 concerning 
the weighted L 1 bound for B±(x,y), the Fourier transforms of m± — 1, which are 
especially needed for low energy estimates. The Marchenco functions B±(x,y) are 
related to m± via 

r"±oo 



/•ItOO 

m±(x,k) = l± B ± (x,y)e 2iky dy. 
Jo 



4.3. Analycity and asymptotics of t(k) and r±(k). The kernel of <fi(H ac ) given 
by (|4.4[) also requires estimates of the coefficients t(k) and r±{k). 

Lemma 4.3. Let V G S. Then t(k), r±(k) € C°° . If \k\ > 1, then 

(4.7) ^ (fe)= / l +^: 1 ) »=° 

v ; dt" w [C^fe-™- 1 ) n > 1 

and 

m 

(4.8) _ r±(fc) = 0(fc — i) n > . 

For the proof of high energy \k\ > 1, we use [TH p. 145] 

1 f 00 

t{ky 1 = 1 ~ 2ik J V ^ m +^ k ) dt 

(4.9) (r±(fc) + l)t(fc)" 1 = 1 =F / h(t,Tk)V(t)m T (t,k)dt. 



For low energy to show £ G C°° the formulas we use are more delicate. 

Lemma 4.4. Let v := VT(0). The following formulas hold, 
a) If v ^ 0, i/ien 
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where W(k) = -2ik + v + V(t)dt J °° B+(t, y)(e 2iky - l)dy. 
b) If v = 0, then 



oo />oo 



t(k)~ 1 = l- V(t)dt (/ B + (t,r])dr))e 2ik Sd£. 

J-oo JO J£ 

Based on Lemma 14.41 we give the proof of the statement t £ C°° as in Lemma 
The statement r £ G°° follows from the fact that t £ C°° and ((43)) . We divide 
the discussions into two cases. 

Case (a) v ^ 0. From the formula in Lemma WM (a) it is easy to see that W(k) is 
C s if J Q y s \B(t, y)\dy < c(t) s+1 , but this is true provided V £ L\ +l according to 
Lemma l4~5l Now since W(0) ^ 0, it follows that t £ C s as long as V £ L\ +1 . 
Case (b) v = 0. From Lemma l4~4l (b) we have if s > 1 

/OO /"OO /-oo 

V{t)dt / / B + (t,rj)drj), 

-oo JO J(, 

which is the Fourier transform of the function 



£^-X(o,oo)(0 / V(t)dt(2iO s B+(t,r,)dr,. 

J-oo J £ 

Observe that in view of Lemma \4. 51 this function is in L 1 if V £ L\ +2 - Hence i _1 
is C s . We conclude that 

m = WF 

is also in C s whenever V £ L\ +2 i since it is a basic fact that v = > 
c > 0, Vfc; cf. e.g., p2 Theorem 1]. □ 

4.4. Fourier transform of modified Jost function. The estimates for m±(x, k), 
t(k) and r±(k), especially in the low energy, depend on the weighted L 1 inequalities 
for the Marchenko functions B±(x,y) and their derivatives. Recall that = {V : 
J(l + \y\V\V(y)\dy < co} and = {V : J(l + \y\V\V^{y)\dy < co,z = 

0,...,n}. 

Lemma 4.5. Let s £ Nq. a) If V £ L^ +1 (IR), then there exists c = c(|| V'|| i i +i ) so 
that for all x £ R 

/oo 
\y\ s \B ± (x,y)\dy < c(l + max(0, T^)) s+1 - 
-OO 

b) If V € W™~ +1 , n > 1, then there exists c — c(|| V|| w ^-i ) so that for all x G R 
\y\ s \d™B ± (x,y)\dy <c(l + max(0,(Tx)) s 

\y\ s \d^B ± (x,y)\dy < c(l + max(0, ( T x)) s . 
The proof exploits careful iterations of the Marchenko equations 

V(t)B + (t,z)dt 



poo p-y poo 

B+(x,y) = / V(t)dt+ dz 

Jx+y JO Jt=x- 

/x+y />0 pt—x+y — z 

V(t)dt+ dz V(t)B_(t,z)dt, 
-oo J y J — oo 
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where we note that B±(x, •) are supported in R±. 

Remark. The inequality in a) is an improvement of [SJ Lemma 3.2, Lemma 3.3], 
where the cases s = 0, 1 were proved via Gronwall's inequality. 

4.5. High and low energy estimates for cf)j(H)(x,y). Applying Lemmas 14.21 
and Lemma 14.31 we can easily prove Theorem 14.11 (a) in high and local energy via 
integration by parts; here the high energy refers to 4>j, j > and local energy refers 
to $ e C °° with supp$ C [-1, 1]. 

For low energy j < 0, (b) of Theorem 14.11 we observe that the condition 
d x m±(x, 0) = for large |x| together with H having no resonance at zero is neces- 
sarily and sufficient for d x <fij(H)(x,y) to satisfy l]4.2[) . Therefore Theorem 14. II (b) 
is true if and only if H has no resonance at zero and m±{x,0) = constant = 1 
for large The latter can occur only when V has compact support in view of 
d x m + 2ikd x m = Vm. 

Remark. Let V E S. If -oo < j < 0, then O still holds for tpj(H)(x,y) with 
£ = 0. However, if H has resonance at zero or V is not compactly supported, 
then (|4.2|) fails for d x (j>j(H)(x,y) in the low energy as we mentioned above. A 
counterexample can be found in [33j for V = —v{y + l)scch 2 .T. For non-smooth 
potentials, in [4] we are able to obtain an appropriate variant of the kernel decay 
(|4.2|) with j e Z, £ = 0, 1 for V = cX[ a ,b]( x ); c > 0, Xe being a characteristic 
function of the set E. 

Problem For non-smooth potentials, it is still open as to determine the class of V 
for which ()4.2|) remains valid especially in the low energy case. 

5. H SATISFYING UPPER GAUSSIAN BOUND 

In this section we prove the analogy of Theorem 11.31 for H with unbounded 
potentials, namely, the Hermite and Laguerre operators. We see that using the 
upper Gaussian bound for Vpt of H we can prove the decay in (|1.2p in a simple 
way, which generalize the results of [HI [21] and [15] . Previously Epperson studied 
the Hermite and Laguerre cases in one dimension using oscillatory integral method; 
later Dziubanski used Hciscnbcrg group technique to prove the kernel decay for the 
Hermite expansion in n-dimension and the Laguerre expansion of integral order in 
ID. 

[55] showed that Assumption 1 1 . 21 is verified when H satisfies the upper Gaussian 
bound (|5.1|) for its heat kernel. The proof is based on a weighted L 1 inequality 
which is a scaling version of [551 Lemma 8]. 

Proposition 5.1. Let I = 0, 1. Suppose V > and H satisfies the upper Gaussian 
bound 

(5.1) \V e x e- tH (x,y)\<c n t- {n+e)/2 e~ clx - yl2/t , Vt > 0. 

If {ipj}j£Z is a dyadic system satisfying (i), (ii), then for each N > 

\W e x ^(H)(x,y)\ <c N y^l 2 {l + yl 2 \x-y\)- N , Vj . 

Remark. The long time gradient estimates for Vpt (t > 1 corresponding to low 
energy) is, in general, not valid for bounded V, not even for positive V € <S(R"). 



HARMONIC ANALYSIS RELATED TO SCHRODINGER OPERATORS 



13 



5.1. Her mite operator H = —A+\x\ 2 . To verifies Assumption 1 1 . 21 it is sufficient 
to show that H satisfies the upper Gaussian bound in (|5.1j) , according to Proposition 

is~n 

For fc € No, let h k be the k th Hermite function with ||/ifc||i,2(R) = 1 such that 

d 2 

(-fa2+ x2 ) h k = (2k + l)h k . 

Then forms a complete orthonormal basis (ONB) in L 2 (M). Let & k (x) := 

h kl (g>---(g> h kn , k = (ki,..., k n ) € Nq- Then {$ fe } is an ONB in L 2 (R"). 
By Mehler's formula Ch.4], the heat kernel has the expression 

e- tH (x,y) = e^ n+2 ^^ k (x)<P k (y) 

fcGNg 

= (27rsinh(2i)) _ "^ 2 e _ ^ coth(2i)(|:c| 2 + |y| 2 )+coscch(2t):r-y 

for all t > 0, x,y £ E". 

One can easily calculated (cf. [55]) that for I = 0, 1 there exist constants c' > 0, 
< c < 1 such that 



\ViPt( X ,y)\<c' 
where pt(x, y) :— e~ tH (x,y). Hence (|5.ip holds 



t -(n+e)/2 e -c\x-y\ 2 /t < t < 1 

e -nt e -c\x-y\ 2 t > i 



5.2. Laguerre operator H = L a = — ^ + x 2 + (a 2 — j)x~ 2 , a > —1/2 in 
L 2 (R + ). The Laguerre is basically a generalization of the Hermite. We will show 
that the heat kernel estimates (|5.ip are valid for L a with a > 1/2. 
For k £ No define the Laguerre function 

where L^(x) are the Laguerre polynomials. Then 

L a M«(x) = (4k + 2a + 2)M£(x) 

and {M^x)}^ is an orthonormal basis in L 2 (M. + ); see e.g. |36] . 

Lemma 5.2. Let a > 1/2. Then the heat kernel of L a satisfies for (. = 0,1 there 
exists constants c'>0,0<c<l such that with n = 1 

(5.2) ^e-^^^^c'l^^,^ ^ l 

Proof. Using the following Mehlcr type formula: if < r < 1, 

oo 

5> fe M^)M, Q (y) =2e-^(x 2/ ) 1 / 2 (l - r)-W 2 e 

fc=0 

x J a (2ixyr 1/2 (l - r)" 1 ) 
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(see e.g. (11) of [21]), we can obtain the heat kernel formula for L a , a> —1/2 

oo 
k=0 

(5.3) =(- J ) Q (sinh2t)" 1 e - £ ^ coth2t (x 2 /) 1 / 2 jy lXV 



V sinh 2i 



where J a is the Bessel function of order a. Now if a > 1/2, (|5.2p follows from the 
usual asymptotics of J a (see [51] or [HI VIII.5]): 



J a (z) 



z a z -> 

— 1 /2 

s ' z — > oo, 



Jz"- 1 z^O 

J aW — S 1/2 



z ' z — > oo. 



□ 

5.3. Laguerre operator in R™ . A natural extension of L a to n-dimensions is for 
a = (ai, . . . , a„), ay > -g 

n 1 

3=1 

in I/ 2 (R™ ) . It can be defined as the Friedrich's extension of the form 

n 

(A,/,/) = £ / l^/^l 2 + + K 2 - j)xj 2 )\f(x)\ 2 dx 

3=1 jR + 

for / G Co°(M"); hence it's domain is a subspace of the Sobolev space VF 2 1 (R rt ). 

The ONB in L 2 («'l) consists of {M^(x) = M^(xi) <g> • • ■ <g> M^{x n ),x = 
(xi, . . . , x n ) G R™, = (fci, . . . , fc„) G Nq }■ Hence the heat kernel of £ Q is given by 

e- t ^{x,y)= E e- tA ^fc(^)^(y) 
fc=(fci,...,fc„) 

2 2 n / ■ \ 

=(-i)l a l( S inh2 i )-"e- M 4^cot h2t -Q (a; . 2/ . ) i/ 2jQj f 

3=1 ^ m ' 

where A£ = X)j=i(4% + 2oj + 2), x, y G R™. Therefore applying the ID result we 
easily obtain the estimates in ()5.2j) for £ Q provided ay > 1/2, j = 1, . . . , n. 

Remark. Thangavelu [47] obtained a heat kernel formula in the special case a 6 No 
via different method. We observe that the formula in [47] Theorem 2.17] for L(a) = 
~d 2 /dx 2 -(2a+l)x~ 1 d/dx + x 2 /A acting in the weighted space (L 2 (R+), x 2a+1 ) is 
equivalent with (|5-3[) for L a acting in the unweighted space L 2 (R + ). 

5.4. Schrodinger operator and associated heat kernel. Proposition l5.1l shows 
that the upper Gaussian bound estimates of e~ tH (x,y) imply the decay estimates 
(|1.2p in Assumption 11.21 However for bounded V ^ the gradient estimates for 
V x pt are not valid in general, this is one reason we work in a more direct way 
to deal with the decay of V x <pj(H)(x,y) for all j G Z as illustrated in the one 
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dimensional case. The heat kernel approach seems to work more efficiently for 
unbounded potentials. For bounded potentials, when the Gaussian bounds are not 
available, we can consider, for instance, the radial case in three dimensions using 
Voltcrra type equation for the eigenfunction of or the non-radial case using 
stationary phase method [361 El 03 EZ] ■ 

6. Conclusion 

The Littlewood-Paley theory of B^' q (H) and Fp ,q (H) surveyed in this paper 
leaves open other problems in the area of harmonic analysis and PDE. Here we 
would like to mention a few problems related to our subject. 

a. The further characterization and identification of Bp' q (H),Fp' q (H), < 
p, q < oo as well as other elements of fundamental theory of function spaces, in- 
cluding distribution theory, have yet to be understood, compared with classical 
Fourier analysis [121 El H2] ■ It may involve semigroup method, Riesz transform as 
well as singular integrals for certain class of rough potentials. 

b. Applying the Littlewood-Paley decomposition to establish dispersion and 
Strichartz estimates for the perturbed wave, Klein-Gordon and Schrodinger equa- 
tions with potentials. Although there has been quite extensive work in this area, 
cf. [31 [71 (37J [5H [53] , regularities of these estimates involving the associated func- 
tion spaces have not received comparable attention. We hope the development of 
the Littlewood-Paley theory in our continuing investigation could give a systematic 
treatment of the regularity problems, which are related to one of the open prob- 
lems concerning Strichartz type estimates for wave equations with potentials in the 
presence of resonance or eigenvalue at zero energy, cf. [3TJ [Ml 131] ■ 

c. For the Laplacian-Beltrami operator — A g on a (complete) Riemannian man- 
ifold (M, g) , in many cases the generic heat kernel estimates are valid (TUJ [TT], [24] . 
However for a Schrodinger operator Hy = — A g + V on M , the heat kernel estimates 
is in general not valid, this makes it unique and more difficult especially for the 
high and low energy analysis of Hy. The treatment here might throw a light on 
considering the analogous problems in a, b on M; in particular we are interested in 
obtaining the analogous results for Hy on Riemannian symmetric spaces. We refer 
to [H [321 EH HQl [45l El [35] for some of the recent development in this direction. 
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